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Abstract 

Let he an N X N random symmetric matrix with independent equidistributed entries 
modulo the symmetry constraint. If the law P of the entries has a finite second moment, it was 
shown by Wigner [T3] that the empirical distribution of the eigenvalues oi Xn, once renormalized 
by -s/iV, converges almost surely and in expectation to the so-called semicircular distribution 
as N goes to infinity. In this paper we study the same question when P is in the domain of 
attraction of an a-stable law. We prove that if we renormalize the eigenvalues by a constant 
Oat of order the corresponding spectral distribution converges in expectation towards a 
law which only depends on a. We characterize and study some of its properties; it 
is a heavy-tailed probability measure which is absolutely continuous with respect to Lebesgue 
measure except possibly on a compact set of capacity zero. 

1 Introduction 

We study the asymptotic behavior of the spectral measure of large random real symmetric matrices 
with independent identically distributed heavy tailed entries. Let {xij, 1 < i < j < oo) be an infinite 
array of i.i.d real variables with common distribution P living in a probability space (fi, P). Denote 
by Xtv the N x N symmetric matrix given by: 

^Nihj) = Xij if i < j, Xji otherwise. 

If the entries have a finite second moment cr^ = IE[a;?j] = / x'^dP{x), and if (Ai,--- ,Xn) are 
the eigenvalues of ^= then Wigner 's theorem (see [14j and generalizations in [9l [1]) asserts that 

the empirical spectral measure YliLi the matrix converges weakly almost surely to the 

semi-circle distribution 



a(dx) = 7tV4:<t'^ — x'^dx. 
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We will consider here the case of heavy tailed entries, when the second moment cr^ is infinite. 
We will assume that the common distribution of the absolute values of the Xij's is in the domain 
of attraction of an a-stable law, for a g]0, 2[, i.e that there exists a slowly varying function L such 
that 

nK\>n) = ^. (1) 
We introduce the normalizing constant ajv by: 

aN = mi{u,F[\xij\>u]<^). (2) 

It is clear that ajy is roughly of order Na ^ indeed there exists another slowly varying function Lq 
such that 

aN = Lo{N)N^. (3) 

We then consider the matrix Aj\f := aJj^X^, its eigenvalues (Ai, • • • , Xn), and its spectral measure 
Aajv := 7! Eili ^\^■ Our main result is 

Theorem 1.1. Let a G]0,2[ and assume ([1]). 

1. There exists a probability measure Ha on M such that the mean spectral measure E,[fiAff] 
converges weakly to /^q,. 

2. ft Am converges weakly in probability to /Xq. More precisely, for any bounded continuous func- 
tion f, j f{x)dfiANi^) converges in probability to j f{x)d^a{x). 

3. Let {Nk)k>i be an increasing sequence of integers such that YlT=i ^ Z*^*" some e < 1, 
then the subsequence fi-Af^^ converges almost surely weakly to fia- 

Remark 1.2. We note that the hypothesis ([T]) concerns only the tail behavior of the distribution of 
the absolute values of the entries. We make no assumption about the skewness of the distribution 
of the entries, i. e about their right or left tails. 

Remark 1.3. It would be useful to control better the fluctuations in Theorem and establish 
almost sure convergence for the whole sequence fiApf ■ 

Our approach is classical. It consists in proving the convergence of the resolvent, i.e of the mean 
of the Stieltjes transform of the spectral measure, by proving tightness and characterizing uniquely 
the possible limit points. We first prove, in section [2l that it is possible , for all later purposes, 
to truncate the large values of the entries at appropriate levels. We then proceed, in section [Sj 
to show tightness for the spectral measures of the truncated and original matrices Ajsf . We then 
introduce, in section HI the following important quantity: for z £ C\M, we define the probability 
measure on C by 

1 ^ 

k=l 

i.e the empirical measure of the diagonal elements of the resolvent of at z G C\M. The 
classical Schur complement formula is our basic linear algebraic tool to study recursively on 
the dimension, as is usual when the resolvent method is used (see e.g [9] or [1]). In section [5l using 
an argument of concentration of measure and borrowing classical techniques from the theory of 
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triangular arrays of i.i.d random variables, we show that the limit points /x^ of satisfy a fixed 
point equation in the space of probability measures on C. Even though we cannot prove uniqueness 
of the solution to this equation, we manage in section [6] to prove the uniqueness of the solution to 
the resulting equation for J x^dfj,^{x), which in turn gives the uniqueness of J xd^^{x). This is 
enough to characterize uniquely the limit points of E[/i^jy] and thus the convergence of E[/i^jy] to 

Ma- 

Once the question of convergence is settled by Theorem ll.H the next question is to describe 
the limiting measure We will discuss in this article three different characterizations of ^a- Our 
approach leads directly to the following first characterization of through its Stieltjes transform, 
defined for z G C\M by: 

G„(z) = j {z - x)-^dfia{x). (4) 
Define the entire function g on C by 

2 f'°° ^ 

ga{y) = - e-^-e-^ydv (5) 



a Jo 

We will also need the constants C{a) = ^"^^ '^('^) ~ cos(^). 

Theorem 1.4. 1. There exists a unique function Y^, analytic on the half plane = G 
C,Imz > 0}, tending to zero at infinity, and such that 

Cia)g{c{a)Y,) = Y,{-zr 

2. The probability measure /x^ of Theorem is uniquely described by its Stieltjes transform 
given, for z £ C"*", by 

G^{z) = -- / e-*e-^(")*^^^dt (6) 



Remark 1.5. Note that fi^ depends continuously on a £ (0,2) since Y^ as described above is 
continuous in a, at list for sufficiently large z, a remark which insures the continuity of Ga{z) at 
list for sufficiently large z and therefore the continuity of fi^ ■ 

Using the characterization given in Theorem 11.41 we prove in section [7] the following properties 
of Ha- 

Theorem 1.6. The probability measure fia of Theorem satisfies 

1. is symmetric. 

2. ^ct has unbounded support. 

3. There exists a (possibly empty) compact subset of the real line Ka of capacity zero, such that 
the measure /Xq, has a smooth density pa on the open complement Uq, = M.\Ka . 

4. Pa has heavy tails. There exists a constant Lq, > such that, when \x\ — > 00 

La 



Pa{x) 
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A second and different characterization of is proposed in the physics hterature by Cizeau- 
Bouchaud [3]. This description has been controversial (see [Ij for a discussion and numerical 
simulations). The strategy used in [3] is also based on the convergence of the resolvent, but on 
the real axis as opposed to our proof of convergence away from the real axis. We unfortunately 
cannot make sense of the strategy used in [3]. We discuss in section [8] the link between our 
characterization given in Theorem 11.41 and the Bouchaud-Cizeau characterization (after correction 
of a small typographical error in [3j already noted by [Ij). 

Remark 1.7. We are not able to prove in general that the exceptional set of Theorem 
is empty, or reduced to zero, even though we conjecture this is true. This would say that /Xq, has 
a smooth density everywhere (except may be at zero) as suggested by numerical simulations and 
accepted by the physics literature. Recent work in progress with A. Dembo indicates that is at 
most the origin for a <1. This question is discussed further in Section^ 



We also describe below (in section [9]) a third characterization of fia, more combinatorial in 
nature. It is based on an extension (due to I.Zakharevich, [15] ) of the classical moment method 
rather than the resolvent approach used both by [3] and us. Obviously because of the heavy tails 
and thus of the absence of moments, one would have to do it first for truncated matrices and then 
try to lift the truncation. More precisely if one truncates the entries at the level Ba^, for a fixed 
B > and define xfj = Xijl\xij\<BaM one can compute the moments of the empirical measures //^s 
of the truncated matrix A^{ij) = a 



N\''JJ ~ "'N ^fj 
k 



X 



df,^s{x) = ^ir{{AU)- 



and study their convergence when N tends to infinity. We establish in section [9] that 

Theorem 1.8. With the above notations, and under the hypothesis of Theorem hl.l\) and the addi- 
tional hypothesis: 

lim ^; , = G 0, 1 7) 

1. E[/i^s ] converges weakly to a probability measure /x^ uniquely determined by its moments and 
independent of the parameter 6. This measure fi^ has unbounded support and is symmetric. 

2. jj^ converges weakly to fi^ as B tends to infinity. 

The moments of /i^ are described combinatorially in Section [9l Thus Theorem 11.81 gives a 
third, independent, description of the limiting measure ^a. As we will see in Section [9l the first 
part of Theorem (jl.Sp is a direct consequence of a general combinatorial result of I.Zakharevich 
and its proof is essentially given in [15]. The convergence of these Zakharevich measures to our 
/id establishes a link between this combinatorial description and the one we have given in terms of 
Stieltjes transforms in Theorem II. 4i This link is far from transparent. 

Remark 1.9. We note that the limiting measure /i^ is in fact independent of the skewness param- 
eter . Thus it is insensitive to the hypothesis about the upper and lower tails of the distribution 
of the entries. This is coherent with Remark \1.2l 
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Remark 1.10. The case a = 2 is covered neither by the classical Wigner theorem (which asks 
for a second moment) nor by our results so far. In fact it is easy to see, using the combinatorial 
approach of Theorem \1.8\ that the limit law is then the semi-circle, even though the normalization 
differs from the usual one. 

Finally, let us mention that the behavior of the edge of the spectrum of heavy tailed matrices 
(when a g]0, 2[) has been established by Soshnikov [I3j. The largest eigenvalues are asymptotically, 
in the scale af^, distributed as a Poisson point process with intensity a~^x~"~^dx. This is in sharp 
contrast with the Airy determinantal process description of top eigenvalues for the case of light 
tailed entries |12j but in perfect agreement with our result about the tail of /i^ given in Theorem 
Ol 



2 Truncating the entries 

Since the entries of our random matrices have very few moments, it will be of importance later 
to truncate them. We introduce the appropriate truncated matrices in this section and show how 
their spectral measure approximate the spectral measure of the original matrices. 

Let us consider (resp. X'^) the Wigner matrix with entries Xijl\;j...\<iBaN for i? > 0, 
respectively Xijl|a,.^.|<^«;a^ for k > 0. Also define 



A 



N 



(Ij^^Xn, a 



A 



N 



a^ ^AT 



Let us remark here that the threshold oat is precisely the scale of the largest entry in a row (or a 
column) of the random matrix X]\[, while the scale of the largest entry (or of the largest eigenvalue) 

9 — 

of the whole matrix is aj^ i.e roughly A^^ . 

We want to state that the spectral measures of the matrices A]\f,Aff and A'^ are very close 
in a well chosen distance, compatible with the weak topology. The standard Dudley distance d is 
defined on P(M) by 



d{n,u) 



sup 

/ll£<l 



fdu - / fdfi 



where the supremum is taken over all Lipschitz functions / on R such that 
norm \\f\\c is defined by 

l/(^)-/(y)l 



£ < 1, where the 



sup ■ 



y\ 



+ sup|/(x)|. 



We will use the following variant di of the Dudley distance. 



di{n,u) 



sup 

ll£<lJT 



fdu - / fdfi 



where the supremum is taken over non-decreasing Lipschitz functions such that ||/||£ < 1 . The 
Dudley distance d is well known to be a metric compatible with the weak topology and the following 
Lemma shows that so is the variant di. 



Lemma 2.1. di is compatible with the weak topology on 
such that there exists /i" G P(R) so that 



, i.e if ^1 is a positive measure on 



lim di{fi'^,n) = 0, 
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then /i" converges weakly to fi and G 'P(M). Reciprocally, if converges to ^ weakly, 
goes to zero. If a sequence fin G V(M) is Cauchy for di, it converges weakly. 

Proof. A compactly supported Lipschitz function / can be written as 

/(x) = /(0)+ r g{y)dy 
Jo 

where 5 is a borelian function bounded by the Lipschitz norm of /. Writing 

fX j-X 

f{x)-f{0)= lg(y)>og{y)dy - \g{y)\lg(^y)<^ody 
Jo Jo 

we see that f can be written as the difference of two non-decreasing Lipschitz functions. Hence, if 
goes to zero as n goes to infinity, / fdfj,n converges to / fdfj, for ah Lipschitz compactly 
supported functions. Hence, fin converges to fi for the vague topology. On the other hand, if /i" 
converges to fi for di, we must have, taking / = 1, 

fi{l) = lim = 1 

which is enough to guarantee also the weak convergence. Indeed, if we now take / G Cb(M), and g 
compactly supported with values in [0, 1], 

|^n(/) - fiif)\ < ll/lloo (^(1) + finil) - fi{g) - fin{g)) + Mg) - fi{fg)\ 

Letting first n going to infinity and then taking g approximating the unit, we obtain the result. The 
second statement is clear since di < d with d the standard Dudley distance (obtained by taking the 
supremum over all Lipschitz functions with norm bounded by one) and the result is well known to 
hold for d. Finally, if a sequence fin is Cauchy for di, it converges for the vague topology (as it is 
tight for the vague topology, and the property of being Cauchy uniquely prescribes the limit) and 
then for the weak topology by the mass property. 

□ 

We next show that truncation does not affect much the spectral measures in the di distance. 

Theorem 2.2. 1. For every e > there exists B[e) < 00 and 6{€,B) > when B > B{e) such 
that, for N large enough 

%{fLA,,fL^B)>e)<e-''^^'^^'^. 

2. For K > 0, and a g]1 — an, 1[, there exists a finite constant C{a, k, a) such that for all N G'N, 

P(di(AA.,/iAx.)>iV^-^)<e-^^"'°s^. 

Remark 2.3. This result depends crucially on the proper choice of the truncation level. Had 
we truncated the entries at a lower level, say N'^a^ with k < 0, then the limit law would be the 
semi-circle. Thus the effect of the heavy tails would have been completely canceled by the truncation. 

Proof. Let X and Y be two N x N Hermitian matrices, and fix and fiy be their spectral measures. 
Then Lidskii's theorem implies (see e.g [8] p. 500) that, if d is the rank of X — y,then 

di{fix,f^Y) < (8) 
Consequently, the following Lemma implies Theorem ()2.2p . □ 
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Lemma 2.4. 1. For every e > 0, there exists B{e) > and 5{e,B) > when B > B{e) such 
that 

¥{rank{XN - X§) > eN) < e'^^^'^)^ 
2. For K > 0, and a g]1 — ok, 1[ there exists a finite constant C{a, k, a) such that for all G N, 

¥{rank{XN - X%) > iV") < e-^^"i°g^ (9) 



Proof, (of Lemma 12.41 ) 

Let Mj~ = 1 (resp. M^^ = 1) if there exists a j < i (resp. j > i such that \xij\ > BaN, and 
Mf = (resp. = 0) otherwise. Define 

N N 

M" = ^ M" and M+ = ^ M+. 

i=l 1=1 

Now let M be the number of non zero rows of the matrix Xjv — X^, obviously 

rank(XAr - X§) < M < M~ + M+, (10) 

so that 

P(rank(XAr - X^) >€N) < F{M~ > ^) + > ^) < 2P(Af- > 

where we observed that is stochastically dominated by M~ (which contains the diagonal 
terms). But if we denote by pi = P(M~ = 1), we have 

Pi = F(3j < i, \xiA > BaN) = !-(!- ^^^"^V < 1 - (1 

where the later inequality holds for c > 1 when N is large enough since 

lim = 1. (11) 

N^oo a% ^ ' 

As a consequence we can estimate the sum 

< AT - ~ NC{B) (12) 

where we denoted An ~ Bn ii An /Bn goes to one as goes to infinity and 

C{B) = l-—{l-e-^). (13) 
c 

For any A > 0, the independence of the M~'s gives 

N N 

E(expAM-) = +pi{e^ - 1)) < exp[(e^ - l)(^Pi)] 
1=1 1=1 
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So that we get the exponential upper bound, for N large enough 



eJV 



IP(M- > — ) < e"-^— E(expAM") < exp[-N<j)_{X,e,B)], 
with 

Obviously, since lim^^oo C{B) = 0, for any e > 0, there exists a B{e) > (of order e""^) such that 
when B > B{e), 

6- (e, B) := sup (A, e,B)>0 

A>0 

and 

eN 

^{M- > — ) < exp[-N5-{e,B)]. 

Using the crude rank estimate (jlOp proves the first claim of Lemma ()2.4p . 

In order to prove the second claim of Lemma (j2.4p . we simply replace B by B{N) = N"^ and e 
by e(A^) = N'^~^ in the proof above. We get then that 

6^{e{N),B{N)) ~ ^(a - l + aK)(iV°-MogiV) 
and similarly for 6-^{e{N), B{N), which proves our second claim. 

□ 

Remark 2.5. We now let A'^ = a^X'^. We note that centering the entries of the matrix A'^ 
defines a perturbation of rank one. Hence, Lidskii's theorem (see shows that 

2 

dl{fLA%,fiA%~E[A%]) < ^• 
Thus we may assume that A'^ is centered without changing its limiting spectral distribution. 

3 Tightness 

We prove in this section that the mean of the spectral measures of the random matrices Ai\f and 
of their truncated versions A^ or A^ are tight. 

Lemma 3.1. 1. The sequence {E,[fiAj^]', N € N) is tight for the weak topology on 'P(M). 

2. For every B > 0, and n > 0, the sequences (E[/1j^b]; N G N) and (ElfiAi^])', N G N) are tight 
for the weak topology on 



Proof. We will use the following classical result about truncated moments (Theorem VIII.9.2 of 
[B]): For any C > a 

/i™ tC-^L{t) C-a ^ ^ 

Therefore, using (|lip . we have 

E[|x.,f lu ,<e,,] ~ (15) 



C-a N 
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or equivalent ly: 



n\A^m'^] - -^B<-"^. (16) 

C, — a iV 



The version for the truncated matrix Afj will also be useful: 

a 



n\Aum'^] ~ 7 — iv'^K-")-^ (17) 

C - a 



Using these estimates with C = 2, one sees that 



sup E[^tr{{A^n ~ (18) 

and that 

sup E[ltr((^^)2)] ~ ^^^^(2-). (19) 

([T2D shows that E[/i^s] belongs to the compact set Kc ■= € 'P(M);^(x^) < C} for any C > 
_2_^2-a g^j^j ]\j large enough. Hence, the sequence (E[/i^s]);A^ G N) is tight, and thus any 
subsequence of E[/i^s] has converging subsequences. We denote by a limit point, i.e the limit 
of a converging subsequence. By a diagonal procedure, we can insure that this subsequence is the 
same for all i? G N, and in particular, since di is compatible with the weak topology, we can find 
an increasing function (p so that for any 6 > 0, Bq < oo, there exists Nq < cxd so that for N > Nq, 
and all B < Bq, 

di(E[/i^s ],/iij) < 5. 

By Lemma [231 and Lidskii's estimate ([8]), we have for all e > 0, 

d,mA,,^,],mA-^J) < E[di(/iA,(,,,A^B^^)] < 26 + e-^(^'^)<^(^) (20) 

with d{e,B) >OiiB> B{e). 

These two inequalities imply that (fis, B £ 'N) is a Cauchy sequence for the modified Dudley 
metric di and thus converges when B tends to oo. Indeed, if we choose e,e\5 > and an integer 
number Bq > B{e) V B{e'), we find that for B, B' G [B{e) V B{e'),Bo] and N > Nq 

di(E[AA,(^J, /iB) <d + 2e + e-'^(^'^)'^W and di{E[fiA,,^^], f^B') < 6 + 2e' + e-^(^''^>(^) (21) 
and therefore 

di{fiB,f^B') < 25 + 2e + 2e' + e-'5(^'^)<^(^) +e-'^(^''^')<^(^). (22) 
Letting going to infinity, and then 6 to zero and Bq to infinity we finally deduce that 

dli^iB,^J■B') <2e + 2e' 

provided that B and B' are greater than -B(e) V B{e'). Hence, fj^B is a Cauchy sequence for di 
and thus converges weakly by Lemma l2. II as B goes to infinity. As a consequence of ()2ip we also 
find that E[/i^_^jj^j] converges to this limit as N goes to infinity. The same holds for the truncated 

versions E[/x^li^^ ]. Thus, we have proved that (E[/i^jy],E[/i^j^])ArgN are tight. □ 
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This lemma can be strengthened into a partial almost-sure tightness result. Consider an 
increasing function : N — > N such that X]Af>o ^{7v) ^ then 

Lemma 3.2. The sequences {fij^B )ArgN, (AA^/jvJAfeN, (Aa**, jA^eN o,re almost surely tight. 
Proof. We note that the truncated moments bound given in (jl6p can be strengthened into a bound 

a 

2" 



in probability as follows. Let M > and C > ^^-B^ Chebychev's inequality reads 



tr((ylf)2) > M + C 



< 




tr((A 



1 ^ 



Ar2 



< 



< 



^rnaxE[<(.,r] 



2aS 



4— a 



a 



M2iV 



(23) 



where we used the independence of the entries at the third step and the truncated moments estimate 
(jl6p for = 4 at the last step. Then Borel Cantelli's lemma implies that for any C > 2^-^^"°^ 

which insures the almost sure tightness of {jj^'^^B^ )N&i- From this point, all the above arguments 

4>(N) 

apply to show the almost sure tightness of {fi'^^^ )NeN and (/i^li^^ )n<^n- D 

4 Induction over the dimension of the matrices 

We borrow the following idea from [3]: in order to prove the vague convergence of (E[/i^^])7veN 
we study the asymptotic behavior, for z a complex number, of the probability measure L^ on C 
given, for / G Cb{C), by 



1 ^ 

k=l 



Here and below, z denotes in short z times the identity in the set of matrices under consideration. 
is thus the empirical measure of the diagonal entries of the resolvent of tIat . In contrast to [3j , we 
will only consider these measures when z G C\M, where everything is well defined since z — Ai\f is 
invertible. 

Note that for z G = {z £ C : 9z > 0}, and for k G {!,••• ,N}, the diagonal term 
((z — AN)~^)kk belongs to the set D := H {x G C : |x| < |9(z)|^-'^}. Lfj is thus a probability 
measure on the compact subset D of C. 
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If we choose the function f{x) = x then 



is the Stieltjes transform of E[/i^j^]. 

Thus, the weak convergence of for all z G C"*" (or even for all z in a set with accumula- 
tion points) would be enough to prove the vague convergence of E[/i^^]. Indeed the latter is a 
consequence of the convergence of its Stieltjes transform, which, as an analytic function on C+, is 
uniquely determined by its values on a set with accumulation points. 

In the following, given a z G C+, we will prove an equation on the limit points of (more 
precisely of its analogue where Aj^i \s replaced by its truncation for some well chosen /c > 0). 
Our main tool will be a recursion on the dimension N, and the Schur complement formula. We 
first investigate how these measures depend on the dimension. 

We let Ajqj^i be the (A^ + 1) x (A^ + 1) matrix obtained by adding to Ajq a first row and a first 
column Aj^{^,k) = ^Ar(/c,0) = a^xok- Hence, ^at+i has the same law as ^at+i- 

We then let Aj^f be the {N + 1) x [N + 1) matrix obtained by adding as first row and column 
the zero vector. 



We also define for z € C\M, 

Gn^i{z) ■= {z - An+i)-' Gn{z) = {z- An)-^ Gn{z) = {z - A^)-^ 

We finally denote by .'^ all quantities where Aj^ has been replaced by its truncated version A 
Thus for z e C\M we define 




'N ~ 



N + 1 



1 




N 



k=l 



N 



'kk 



k=l 



Lemma 4.1. 



1. G'j^{z)kk is equal to G'j^{z)kk for k > I and to z ^ for k = 0. 



2. 



hm ^E[|G^+i(^),,-G^(z)fc,|]=0. 



k=\ 



3. For K g]0, 2r^[ and < rj < i(l - k{2 - a)), 



lim P(d(L^^L^';,)>7V-'') = 



Here, as above, d is the Dudley distance on V{C). 



Proof. We note that 




(24) 
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which immediately yields the first point. For the second, let us write 

AT 



kk 



1=0 



N 



+ X] G'Af+l(^)fco(^Ar+i - AN)oiG%{z)ik 
l=Q 

N 

G^^+^{z)koY.A%mG'h{z)ik 



1=0 



where we noticed above that G'j^{z)Qk is null for /c 7^ by (j24p . Therefore, we find that 

N 

nG%+i{z)kk-G'h{z)kk\? < E[\G%^,{z)ko\^]E[\Y,A%mG%{z)ik\'] 

1=0 

by Cauchy-Schwartz's inequality. We recall that we have seen in remark 12.51 that we can assume 
that the entries of the matrix A'^ are centered. Using then the independence of Aqj^ and Gm{z), 
summing over k G {1, • • • , A^} and with a further use of Cauchy-Schwartz's inequality, we find that, 

N 



^Y.K[\G-^^^{z)kk-G%{z)kkW 



k=l 



< 



maxEp^(0i))2]l -J]E[|G^^_i 



Z)kQ\ 



E 



k=l 



1 ^ 



N\^)lk\ 



Lk=l 



We now note that the entries of the resolvent Gn{z) are uniformly bounded in modulus. Indeed 
observe that, if [/ is a basis of eigenvectors of A^, with associated eigenvalues (Aj, 1 < i < A^) G M^, 
for any A;, / G {O,--- ,NY, 



\G%{z) 



kl\ 



r 



9(z)| 



(25) 



and the same holds for Moreover, since the spectral radius of Gn{z) is bounded above 

by 1/1 9(2;) I, we also have 



1 ^ 1 

- \GUz)ik\' = j^t<GUz)G%{zr) < 



N + 1 



iV|?>(z)|2' 



Lk=0 
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Hence, we deduce 



N 



k=l 

But we know how to control the truncated moments E[(74^(0Z))^]. Indeed by the estimate (jl7p we 
see that there exists a finite constant c such that for all N £N, 

nrax EP^(ij)|2]<ciV- (26) 

l<i<j<N 

with e = 1 — k{2 — a) > 0. The proof of the second point is complete. 
We finally deduce the last result simply by 

1 ^ 

k=0 



1 N + 1 1 . 



and since Gn{z) and Gm{z) differ at most by a rank one perturbation, 



... 1 , 

1 



= FTT^ + FT! 



implies that 

2 



d{L^jf,L%'') < 



N + 1 

This shows by Chebychev's inequality that for all ?] < | 



N^oa 



□ 



To derive an equation for L^j^'^, our tool will be the Schur complement formula, which we now 
recall. Let ^at+i and be as above. 

Lemma 4.2. For any z S C, 

-1 

((^;v+i - ^^)"')oo = I ^^(00) - ^ - E ^NmANilO) {{An - zl)-^) 

\ k,l=l 

Proof. The proof is a direct consequence of Cramer's inversion formula: 

. - 1. det(^jv - zIn~i) 

{{An^,-zI) )oo= det(i^+i-z/) • 
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To get a more explicit formula for this ratio, write 



A 



N+l 



zl 



^Ar(OO)-^ Oo 



An - z 



with Oo = (A(01), • • • , A{{)N))^ and use the representation 





' I - 


-BD-^ ■ 




' A 


B ■ 









/ 




C 


D 




z, B 


= ao 


C = al 


and D 


= A 


N — 



A-BD-^C 
C D 



we conclude that 

<lei{AN+i - zl) = dei{AN - zI)(let[A{m) - z - {ao,{AN - zI)~'^ao)] 
This proves the lemma. 



□ 



We now show that, in the Schur complement formula above, the off-diagonal terms in the sum 
in the right hand side are negligible. 

Lemma 4.3. For any 5 > {), any z with \'^{z)\ > 5, any < «; < 2{2-a) ' ^"^^ R> 

2 



Y^AUOk)AUOl){{A%-zr')J>R) 

k^l 



< 



^2^2.-1^2- 



with e = 1 — k{2 — a) > ^. 

Proof. Following Remark 12.51 we can always assume that the entries of are centered. By 
independence of A'^{Qk) and A^, we find that the first moment of the off-diagonal term vanishes: 

E[ ^ A%{Qk)A%m{{A%-zI)-^)J = Q 

l<k,l<N 



and that the second moment is small: 



'kl 



< 



2(maxE[(^^(ii))2])2E[^ | ((A^ - z)-^)^^ p] 



kl 



< 2iV-^^E[tr((A^ - z)-\A% - z)-')] < 



c>/'^^|2 • 



Chebychev's inequality concludes the proof. 



□ 



We finally derive from the previous considerations a first approximation result for L^^^ . This 
will be our first step to obtain a closed equation for the limit points of the spectral measure (such 
an equation will be derived in the next section). 

Lemma 4.4. For < k < 2{2~a) ' e = 1 — k{2 — a) > ^.Let z G C^. For any bounded Lipschitz 
function f , 



lim \nL'^\f)]-¥. 

Ai— >oo 



f\{z-Y,A%m?GUz)ki)j 



-IN 



0. 
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Proof. It is clear, by Lemma 14.11 that it is sufficient to prove that, for a constant c',and every 
Lipschitz function / 



|E[L^;,(/)]-E 



N ^ 

f\ \z-Y,A'k{OkfGUz)kk 

k=l / 



< 



c\\I\\c 



\Q{z)\iN— 



(28) 



We have proved above that, for z G C\M, there exists a random variable e]\f{z), the sum of the 
off diagonal terms and An (00) 



\\eNiz)\ >R)< 



+ 



4a 



^2jv2e-l|cj(2)|2 R2(2-a)N 



such that 



N 



^^+1(2)00 = [z-Yl ^NiOk)^G%{z)kk + eN{z) 

V k=l J 

In particular we have for any Lipschitz function /, 



N 



/I \z-Y,^Nm''G%{z\u^en{z) 

k=\ J 



E[/(G^+i(^)oo)]=IE 
Observe that with A^^ = C/diag(A)C/*, 

G%{z)kk =^\uki\'^{z- Xi)-^ 



N 



1=1 



is such that 



Q{z)Q {Gg{z)kk) < 0, \G%iz)kk\ < mz)\-^. 
In particular, we always have 

^(^-Ek=iAf,m'G-^{^h, 



Qiz) 

Thus, on |eAr(z)| < |9(2;)|/2, we obtain the control 

z - A%iOkfG%{z)kk + SNiz)^ - (^z 
Hence, if / is Lipschitz, 



> 1. 



N 



-1 



Y,A%m'G%{ 



Z)kk 



k=l 



(29) 



< 



2|£jv(z)| 
|9(.)|2 • 



N \ 



E[/(G^+i(^)oo)] =E[f{[z-Y, AUOkfGUz)kk ]+0(ll/lk)(IE[^y Al] + 



k=l 



e N{z)\ 



where the last error comes from the weight of putting and removing the cutoff |eAr(z)| < |9(2;)|/2, 
due to the fact that ||/||£ also bounds the uniform bound on /. Now, the right hand side does not 
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depend on the choice of the indices and so we have the same estimate for all K[f{G']^^^-^{z)i.k)], for 
k € {0, 1, • • • , A^}. Summing the resulting equalities we find that 



IE[i^+i(/)]=IE 



N N 
k=l / 



1 



X5Z 



Af2.-l|9(^)|4- 



This proves the estimate (I28p and thus the lemma. 



□ 



5 The limiting equation 

We prove in this section that the limit points of the sequence of measures E[L^'^] satisfy an implicit 
equation. This section will rely heavily on a result about the convergence of sums of triangular 
arrays to complex stable laws. We have deferred to Appendix [TU] the statements and proofs of 
these convergence results. We also refer to the same Appendix for notations and references about 
complex stable laws. 

Hereafter z G will be fixed. We have seen that IE[L^'^] is a compactly supported probability 
measure on C (since its support lies in the open ball with radius 1/|9(2)|). Therefore, (E[L^'''])ArgN 
is tight, and we denote by //^ a limit point. Recall that for z G C"*", /i^ is a probability measure on 

c-n{\y\<i/mz)\}. 

In order to state the main result of this section we will need the following notations. For t, z £ C, 
we denote by {t, z) the scalar product of i and z seen as vectors in R^, i.e {t, z) = 
For a probability measure on C, and t G C, we define the numbers cT^,Q(t) and (3^,a{t) by: 

a,,o.{t) = [^J\{t,z)rd^,{z)]^ (30) 

and 

_ ! \ <t,z > I'^sign <t,z> d^jz) 
j \ < t,z > I a/i(2) 



where 



c'-i fs^.ittm ,3.) 



Definition 5.1. For a probability measure jj. on C, we define the probability measure P^' on C by 
its Fourier transform 



I 



e'<''^>dP^ix) = exp[-a^ .(t)t(l - .(t)tan(^))] 



Pf^ is well defined by this Fourier transform, indeed is a complex stable distribution. For 
this description of P^ see the appendix [TOl 

We can now state the main result of this section. 

Theorem 5.2. For < k < 2(T^^> limit points fi^ o/E[L^'^] satisfy the equation 

I fd,,^ = j /(j^) dP^^\x) 
for every bounded continuous function f. 
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Proof. We consider a subsequence of (E[L^'^]) converging to fi^, i.e an increasing function (j){N) 
such that (E[L^'|-'^^]) converges weakly to fi'. We denote by P^ the law of J2k=i(.ANiOk)fG%{z)kk- 
For z G C"^, P^ is a probability measure on C~ since then GN{z)kk G for all k. If f is Lispchitz, 
Theorem 15.21 is a direct consequence of the main result of the preceding section, i.e Lemma \AA\ 
and of the next crucial Lemma [ 



Lemma 5.3. /f E[L^'j.'jy^] converges weakly to fj,^ as N goes to infinity, then -P^(Ar) converges weakly 
to P^' as N goes to infinity. 

It is then easy to see that the statement of Theorem 15.21 extends to any bounded continuous 
function. □ 



We now have to prove Lemma (j5.3p . 
Proof. We apply first the following concentration result for i^'^- 

Lemma 5.4. For k £ (0, 2^)' let e = I — k{2 — a) > 0. There exists a finite constant c so that 
for z £ C\M and any Lispchitz function f on C 



9(z 



This Lemma shows that since E[L^'|.'jy^] converges weakly to n^, then -^^(•'^■j also converges almost 
surely to the non random probability /x^. From there, one can apply Theorem 110.31 of Appendix 
[TUlor more precisely its extension Theorem 110.41 which has been built to fit exactly our needs here, 
when applied to the variables = A{0, k)"^ and therefore a^ = a^. One must simply notice that 
the exponent a in Theorem 110.41 must be replaced here by ^. This concludes the proof of Lemma 



Proof of Lemma ( (5.^[ j . We prove this concentration lemma using standard martingale decomposi- 
tion. We assume that / is continuously differentiable, the generalization to any Lipschitz function 
being deduced by density. We put 

1 ^ 

FM{A%{kl), k<l):= L%''{f) = ^ E fiGNiz)kk) 

k=l 

Let n = N{N -l)/2 + N and index the set {A%{kl), k < I) hy {Af, l<i< N{N -l)/2 + N for 
some lexicographic order. Then, if we let J^i = a(Aj,l < j < i), the independence and identical 
distribution of the Af^s shows that, if Pn denotes the law of A^ (i.e the properly truncated and 
normalized version of P), 
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E[(F^-E[Fjv])2] 



n-l 



Y,mnFN\J'^+l]-E[FN\H? 

i=0 



i=0 
n-l 



< ^ j {Fn{xi,- ■ ■ ,Xi+i,- ■ ■ ,Xn) - j Fn{xi,- ■ ■ ,Xi,y,Xi+2- ■ ■ Xn)dPN{y)fdP^"-{x) 

1=0 

n—1 „ 

< Y^W^-'^^^^^Wlo U^-yfdPf{x,y) (33) 

1=0 

In our case, for all A: G {1, • • • , N}, all m,l £ {!,■ ■ ■ , N}, 

dA^J{GN{z)kk) = f {GN{z)kk){GN{z)klGN{z)mk + GN{z)krnGN{z)lk) 

which yields 



1 ^ 

dA^lFN{A) = —^f{GN{z)kk){GN{z)klGN{z)mk + GN{z)kraGN{z)lk) 



k=l 

= ^ {[GN{z)D{f)GN{z)U + [GN{z)D{f)GN{z)U) 

with D{f') the diagonal matrix with entries {f'{G]\f{z)kk))i<k<N- Note that the spectral radius of 
GN{z)D{f')GN{z) is bounded by ||/' ||oo/|9(^)P and so sinc~e for all /, m e {1, • • • ,N}^ 

\[GNiz)Dif)GNiz)U\ < ||G;v(^)I)(/')G'iv(^)||oo < l|/'||oo/|9(^)|' 
we conclude that for all l,m £ {1, • • • , A^}^, 

Thus, ([33]) shows that 

^ ^11/ lloo ly-t 



where we used the truncated moment estimate (117p . Chebychev's inequality then provides the 
announced bound. 

□ 
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We now apply Theorem 15.21 for a particular choice of the function f. To this end, we need to 
define, for any a > 0, the usual branch of the power function x — > x", which is the analytic function 
on C\M~ such that (i)" = . This amounts to choosing, if x = re*^ with g] — 7r,7r[, 

This function is analytic on C\M~ and extends by continuity to a; = re*^ with 9 decreasing to — tt 

lim (re'^T = r^e"^™ 

When X = re*^ is on the other side of the cut M", i.e when is slightly smaller than vr, the function 
jumps by a multiplicative factor e^*"'^. We want to choose in (|5.2|) the analytic function /(x) = x^ . 

Theorem 5.5. ForO < k < 2{'2-a) ' ^ limit point o/E[L^'^] and define X^z := J x^dp,^{x). 

Then 

1. X,,z is analytic in and \Xz\ < — ^— o- 

2. Xfj_z is a solution of the following equation: 

l-OO 

X^.=iC{a) (ft)t-^e**^exp{-c(a)(it)tx^4dt. (34) 
Jo 

with C{a) = and c(a) = r(l — 
Proof. The first point is obvious. Indeed, for some increasing function (j), 

k=l 

For each A^, Xf^ is an analytic function on C^. Moreover, \X^\ < — ^—a for all N. This entails 

'2 , [ z [ - 

that any limit point X^z must also be analytic in C^. 

In order to prove the second point and obtain the closed equation (j34p we will need the following 
classical identity: 

Lemma 5.6. For all z G C"*", 




/•oo 

= iC{a) / {it)'^-^e'^'dt 
Jo 



with C{a) = 

This Lemma is proven by a simple contour integration, it is also a consequence of Lemma 16.21 . 
proven in the next section ( plug y=0 in the statement of Lemma 16. 2p . 

By Theorem (|5.2p . and since /x^ and are supported in C~, we can write 
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Applying Lemma 15.61 to z ^ z — x ^ for P^^" almost all x, and integrating over the x's we 
have, by Fubini's theorem, 

X^.=iC{a)j {it)^-^e'^^ j e-^^''dP^'\x)dt. (35) 

We now use Theorem 110.51 in the appendix, with v = ^i^ here, and replacing a in Theorem 1 1 . 5 1 bv 
f . We see that: 

j e-^*^dP^'(x) = exp{-c(a)(it)t j x^dfi'{x)}. (36) 
Plugging this equality into ([351) yields 

X^^=iC{a)j {it)'^-^e'^'exp{-c{a){it)'^ J x'^dfi'{x)}dt. (37) 

We have obtained the announced closed equation 

POO 

X^.=iC{a) (ii)t-^e**^exp{-c(a)(it)tx^4(it. (38) 
Jo 



□ 



6 Proofs of Theorem 11.11 and of Theorem 11.41 

In this section we gather the preceding arguments and prove Theorem 11.11 and Theorem II. 4i This 
proof will be based on the following uniqueness result for the closed equation (jMj) . We recall the 
notation 

2 f'°° ^ f'°° 
ga{y)-=- e-^'^e-^ydv = d-h'^exp{-dy}dt 
a Jo Jo 

Theorem 6.1. 1. There exists a unique analytic function of z & , such that \Xz\ = 
0(|/m(z)|~ 2") at infinity, satisfying the equation 

l-OO 

X, = iC{a) I (it)t-ie**^exp{-c(a)(it)txjdt. (39) 
Jo 

2. This solution in fact also satisfies: \Xz\ = 0(|2;|~2"). 

3. If one defines Yz := {—j)^Xz, then is the unique solution of the equation 

{-zrYz = C{a)g^{c{a)Yz). 

analytic on C+ and tending to zero at infinity. In fact =0(|z|~") 

Proof. We already know that there exists such an analytic solution X^. Indeed we have seen in 
the preceding section that, if /^^ is a limit point, then X^z is such a solution. In order to prove 
uniqueness, we will use that: 
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Lemma 6.2. For all z G C"*", and any y G 

1 



oo 



(— = i / (it)^"^e*'^exp[-(-z)2(it)2y]dt 

^ JO 

Proof. 

We write 2; = re*^ with some G]0, 7r[. Assume first that e]0, |[. Since /(n) = (u) 2'~^e"^e~"^ 
is analytic in C\M^, for all i? > finite, its integral over the contour 

r = {it, e<t<R}U {e'"R, 77 G [|, vr - 0]} U {e'^~'h, R < t < e} U {e'\, r, e [tt - 9, ^]} 

vanishes. Note that r] + 6 £ + 9,-k] so that 3f?(i?e'''z) = Rrcos^r] + 9) < for ah rj £ [f ,7r - 6*] 
and 9 e]0, §[. 

This shows that 

lim Rf{e"iR) = V?? G [-, vr - 0] ^ lim / f{e'''R)dr] = 0. 

Similarly, 

lim sup I / /(e*''e)(i7?| < cx) ^ lime / /(e*''e)dr/ = 

Hence, letting — > 00 and e — > 0, we find 

i / f{it)dt + / /(e^("-^)t)e^("-^)dt = 0. 

JO J+00 

In other words. 



where we finally did the change of variable t' = \z\t. 

Noting the obvious facts {—z~^t)^~'^ = {—z~^)^~^t2~^ and {—z~^)^^^{—z)^~^ = 1, we thus 
have proved that 

/.oo poo a 

i / (it)t-ie^*"exp{-(-z)t(ft)ty}dt = (-z-^)t / -le-VJ'^^dt 
Jo Jo 

which proves the claim when 9 g]0, The case 9 £ [^,7r[ is identical after an immediate mod- 
ification of the definition of the contour. It can also be obtained by a trivial analytic extension 
argumentn 

By Lemma 16.21 we remark that, if is a solution of the equation (j39p and if z = \z\e^'^, 

/•oo 

X, = -e-'^C{a) / {-e-'^t)'^-^e-'\'\ exp{-c{a){e-'H)'^ X,}dt (40) 
Jo 

= —C{a) / (— )f-ie~*exp{-c(a)(--)fX,}dt 
z Jo z z 

= ( — )tC7(a) / tt-ie-*exp{-c(a)tt(-i)fx,}dt. (41) 
2 Jo ^ 
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Hence, if Yz := (— ^) 2 Xz, we obtain 

POO 

{-z)'^Yz = C{a) / tt-ie-*exp{-c(a)ttyj(it. (42) 
^0 



This equation for can be written simply as 

{-zrYz = C{a)g{c{a)Yz). 

We recall that we have assumed that there exists a constant Ci such that \Xz\ < Ci9(z)~ 2". 

Now, consider the function of two complex variables F{u, y) = uga{y)—y. Obviously -F(0, 0) = 
and 5y-F(0, 0) = —1. By the local implicit function theorem, there exists ei > and £2 > 0, such 
that for every u G C with |m| < ei there exists a unique y{u) G C with \y{u)\ < €2 satisfying the 
equation F{u,y{u)) = 0, i.e uga{y{u)) = y{u). Moreover 

\y{u)\<C\u\. (43) 

For any z £ C+, such that 9(z) > L, with > ^^77^, tlien \Xz\ < CiL~t so that 

\Yz\ < ji < j^. Thus for z G C+, such that Q{z) > L we have that 



' < ei, \c{a)Yz\ < €2 



'C{a){-z 

Thus the uniqueness in the local implicit function theorem shows that Y^ is given by 1^ = 
c(ay^( c(a)(-z)° ) ^^'^ ^^^^ ~ ( — iyi"^^' ^^^^^ "'^z analytic on 2; G C"*" and uniquely 

determined on the set of 2; G C'^ such that 9(-z) > L it is uniquely determined. This proves the 
claim of uniqueness for Xz- Using the bound (j43p now proves the improved bound at infinity, i.e 
\^z\ = 0(|z|~ 2"). These arguments prove the second and third statements of the theorem. □ 

We can now deduce from this last uniqueness result the convergence of the mean of the nor- 
malized trace of the resolvent. 

Theorem 6.3. For any k g]0, 2(2-0) i> ^"^V ^ ^ > ^[Tf'^k=i^N(^)kk] converges as N goes to 
infinity to 

G^{z) := i r e**^e-^(")(^*)*^Mt = -- H e-'e-^'^^'^^^dt (44) 



z Jo 

Proof. For any z G C+ and any limit point ^u^ , 
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The uniqueness of implies that the mean of the resolvent E[A^ ^ii{z — A'^) ^] has a unique 
limit point which is given by 

Jo 

This shows that E[N-hi{z - A%)-^] converges to Gaiz). In order to finish the proof, observe that 
for z E C"*", we can use the same arguments than in the proof of Lemma 16.21 to see that 

G„(z) = i / e'''e-<''^^^^dt 
Jo 

z Jo 

= _i /"e-'e-^W^^^di (45) 
z Jo 



□ 

This last result enables us to conclude the proof of Theorem 11.11 and Theorem II. 4[ 
Proof of Theorem and Theorem \1.4\ 

By Lemma Em E[/i^7^] is tight for the weak topology. Taking any subsequence, we see that 
any limit point /i is such that its Stieltjes transform must be equal to Ga{z) for all z G C^. This 
prescribes uniquely the limit point /i and thus insures the convergence of E[/i^7^] towards /i G 'P(M) 
so that 

j\z- x)-^dfi{x) = Ga{z), z e C+. 

By Corollarv 12.21 and the fact that 

di(E[/iA^],E[AA^]) < E[di(/iA^, Aa^)] 

we also conclude that E[/i^j^] converges weakly towards /i. By Lemma 15.41 for any z E C\M, 
L^'^(x) = J{z — x)~^dflA'^{x) converges in probability towards Ga{z). This convergence holds as 
well for finite dimensional vectors (J {zi — x)~'^dfiA'^{x), 1 < « < n). Since {{z - e C\M} is 

dense in the set Co(M) of functions on M going to zero at infinity, we conclude that J f{x)dfiAi^{x) 
converges in probability towards / f{x)d^{x) for all / G Co(M). But also fiA'^ (1) = = 1 and so 
this vague convergence can be strengthened in a weak convergence (see the proof of Lemma l2.ip . 
We finally can remove the truncation by k by using Corollarv 12.21 Again by Lemma [5.41 L^''(x) = 
J{z — x)~^dfiA;^{x) converges almost surely along subsequences (/)(iV) so that ^(/>(A^)~'^ < oo 
by Borel-Cantelli Lemma. As e = ^ — is as close to one as wished, for any sequence (l){N) 
so that ^(/)(A^)~^ < oo for some e < 1, we can choose k close enough to one so that L^'^'jy^(x) 

converges almost surely to G{z). This entails the almost sure weak convergence of A^j^^) 

same arguments as above. □ 

Remark 6.4. // we could prove that the equation given in Theorem \5.2\ admits a unique solution 
fj,^ , at least for z in a set large enough, the convergence o/E[L^'^] to this solution would be assured. 
We cannot prove this uniqueness result. But as we have seen we do not really need such a strong 
uniqueness statement either. We rather have proved a weaker statement, i.e the uniqueness of 
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/ xd^^{x), which already entails the uniqueness of the limit points for "Klj xdL^j^{x)], i.e the mean 
Stieltjes transform of the spectral measure of A'^. This is sufficient for our needs hut the question 
of the uniqueness of solutions to the equation given in Theorem \5.^) remains intriguing. 

7 Study of the limiting measure. Proof of Theorem 11.61 

In this section, we prove Theorem 11.61 First, the fact that the Hmit measure Ha is symmetric is 
obvious. It suffices to consider the case where the entries have symmetric distributions. To prove 
the other statements, we need to consider the Hmit of Ga{z) as z tends to a positive real number x. 
We first remark that the analytic function defined on C"^ is univalent (i.e one-to-one). Indeed 
this is an obvious consequence of the equation, valid for z G C^: 

{-zTY, = C{a)go,{c{a)Y-,). 

In order to study the boundary behavior of Gct{z), we thus have to study the boundary behavior 
of the univalent function Y^. For x G M, the cluster set Cl{x) is defined as the set of limit points 
of Yz when z tends to x (see [5] or [10]). It is easy to see that for any non zero x G M the cluster 
set Cl{x) is reduced to one point in C U {oo}. Indeed, assuming w.l.o.g that a; > we have, for 
any finite v G Cl{x), the equality C{a)gaic{a)v) = e*™(x)°f . If Cl{x) contains two points it is a 
continuum, i.e a compact connected set with more that one point (see [5]). By analytic continuation 
we would then get the equality C{a)ga{c{a)v) = e^'^"{x)"v for every v € C which is false. The only 
remaining possibility for Cl{x) is to be reduced to one finite point or to the point at infinity. We 
define 

= {x G R,C/(x) = {oo}] 

We first prove that K'^ is bounded. The proof of Theorem 16.11 using the local implicit function 
theorem at infinity, shows that Y^ admits and analytic extension to the setjz G C, \z\ > L} for L 
large enough, and that this extension satisfies = 0(|2;|~"). This obviously proves that, when 
|x| > L, the cluster set Cl{x) is reduced to one finite point and thus that K'^ is bounded. 

We consider the complement f/^ of . Let x £ U'^ and Yx the unique point in the cluster set 
Cl{x). By continuity, for x > 0, satisfies the equation 

e-^™x-y, = C{a)g{c{a)Yx). 

The local implicit function theorem can be applied to this equation at (x,!^), except for the 
subset say F of M where the derivative vanishes. The exceptional set F must be bounded, since the 
derivative does not vanish at infinity, and its points must all be isolated. Thus F is finite. For any 
X G U^\F, the implicit function theorem shows that Yz can be extended analytically on a complex 
neighborhood of x. Hence Ua '■= U^\F is open and its complement Ka = K'^ U F U {0} is closed. 
Ka is also bounded and thus compact. 

Finally we use Beurling's Theorem which states that the set K'^ has capacity zero, and thus 
also the set Ka (see [5] or pO]). 

For any point x in the open set Ua the function Yz admits an analytic extension to a complex 
neighborhood of x, and thus the Stieltjes transform Ga{z) admits a smooth extension, which proves 
that ^a has a smooth density pa on the open set Ua- Indeed, for x G Ua 

limGa{z)= Hfia (x) - inpa (x) = / e'^e-^^")* ' dt 

z^x X Jq 
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In particular the density of the measure fia is given, if 1^ = r^e^'^'' , by 

p^(x) = — e-*e-"(°)*^["-'=°'('^-)lsin[c(a)ttr^sm(</.^)](it. (46) 



TTX 







Note that we now know that is well defined and smooth for x large enough. We also have seen 
that Yx = 0(|x|~") and thus that Y^ ~ e*™C(a)(7a(0)x~" . Hence, when x — > oo, the following 
asymptotic behavior holds for Ga{x) = liniz^x Gaiz) ■ 



X 



1 1 /"OO 

-/ e-^{l-ct^Yx{l + o{l)))dt^-{l-c e'H^ dtYx{l + o{l))) 
X Jo X Jq 

Identifying the imaginary parts of both sides we get: 

Pa(a;) ~ ci (a) . 

X 

Which proves the last statement of Theorem 11.61 

8 Cizeau and Bouchaud's characterization 

In [3], the authors propose the following argument; they look at Gn{z)qq for z on the real line. By 
arguments similar to those we used (but with no a priori bounds on the Giy{z)kk) they argue that 
G]sr{z)oQ converges in law as N goes to infinity. The limit law, that we will denote Pq to follow 
their notations (but which is /x^ in ours) is then given by the implicit equation (11) in [3j 

/ f{y)dPG{y)= [ f{^)dPs{y)= [ ^f{y)dPs{z--). 
J J z-y J y 

Ps = La^^^'^^^'^ is now a real-valued stable law with parameters G{z) and P{z) given self-consistently 
(see (12a) and (12b) in [3]) by 

C{z) = l\y\'^dPG{y) = l\y\^-^dPs{x-^) 

P{z) = / \y\^sign{y)dPG{y) 



where there was a typographical error in the definition of /3 in 12b which was already noticed 
in [3]. We in fact have that for any real t, 

e-'^^dPsiy) = e"'^^'**('^('=)"^*'"'('^)'^('^)) 



where we used that := / (x) 2 dPcix) = e 4 [cos(^)C(2;) — ism{^)P{z)]. So, we see that the 
description of the limit law is very similar to ours, except that z is supposed to belong to M. Let 
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us assume (as seems to be the case in [3]) that C{z) and P{z) are finite. Then, also is finite 
and we see that for non negative real z's 

= Jiz-yr^dPsiy) 

= -C{a) / tt-ie-*^e-r("~i)(^*)^^^dt. (48) 



Hence, Kz and the Xz introduced in section [7] satisfy formally the same equation, except that 
Xz satisfies it for z G C"*" and Kz for real z's. Moreover, we have seen that Xz can be extended 
continuously to z real in {K'^Y and then this extension Xz satisfies the same equation that Kz- 
This indicates that we expect Kz and Xz to be equal, at least on {K'^y. In fact, Xz is the unique 
solution of this equation with an analytic extension to C"*" and going to zero at infinity. In [3], 
under (12a-12b), it is claimed that the equations defining C{z),l3{z) have a unique solution, and so 
that Kz is also determined uniquely by (08]) . We could not prove the uniqueness of the solutions to 
this equation on the real line. In any case, if we beleive either that Kz extends analytically on C"*" 
and goes to zero at infinity or that the above equation has a unique solution for z G M, we must 
have Xz = Kz at least for z G {K'^y. 

The second claim of |3] is that the density of the limiting spectral measure p{z)dz = dfi{z) is 
given, see [3] (14), by 

/ X dPs f X 
dz 

Note that by Fourier inversion, if Kz = Xz, for z > 0, since Ps is a probability measure on M with 
Fourier transform given by ()47p . 



dz 2tt 
1 



oo a 



e 



IT \ Z 



oo 







and therefore we miraculously recover our result ()46p . Hence, at least for z G (K'^Y, the prediction 
of [3] coincides with our result if we beleive that (jlHj) has a unique solution. 

9 The moment method. Proof of Theorem 11.81 

We prove here Theorem 11.81 using the moment method developed by I. Zakharevich [15j. For any 
> 0, we consider the matrix X^ with truncated entries = Xijl\ri;-j\<Bapf and the normalized 
matrix Af^ = a~j^X^. Recall that work here under the additional hypothesis ([7]): 

lim^^44^ = ^G[0,l] 

We begin by the following estimate on moments of the entries of A^. 
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Lemma 9.1. For any integer m > 1, the following limit exists 



Cm. = lim 



E[A^{ijy 



Moreover, if m = 2k is even 
If m = 2k — 1 is odd 



N-.00 iV--iE[yl^(ij)2]- 



m — a a 



m — a a 



Proof. It is a simple application of the classical result about truncated moments (Theorem VIII. 9. 2 
of ^ ) already used in Section 13.11 (fT5l) : For any ( > a 

The first item of the Lemma is a direct consequence of this estimate for C = 2 and C, = 2k. The 
second is also a consequence of this estimate, used for x{ij)~^ and x{ij)~ , and of the additional 
skewness hypothesis ([7]). 

□ 

This lemma enables us to get the main result of this section, i.e the convergence of the moments 
of the spectral measure of the matrix A^. We will need some more notations that we take verbatim 
from Zakharevich. For any integer A; > 1, we define as the set of all (ei, ei)such that Yli=i — 
k and ei > €2 < ... > e; > 0. For any (ei, ei) G Vk define T{ei, e/) as the number of colored 
rooted trees with k + 1 vertices and I + 1 distinct colors, say (ci,...q) satisfying the following 
conditions: 

1. There are exactly nodes of color Cj. The root node is the only node colored cq 

2. If nodes a and b are the same color then the distance from a to the root is the same as the 
distance from b to the root 

3. If nodes a and b have the same color then their parents also have the same color 

With these notations we have the following convergence result, directly implied by Zakharevich's 
results. 

Lemma 9.2. 1. For every integer k > 1, the following limit exists 

hm E[/ x'dftMx)]=:m^ (49) 

2. m^ = Oifk is odd, and mg, = E(ei,...,eOeli. ^(^i^ TlLi C2e, • 

3. There exists a probability measure fi^ uniquely determined by its moments . fi^ is inde- 
pendent of the skewness parameter 9. 

4. Ha has unbounded support and is symmetric. 
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5. The mean spectral measure E[/i^s] converges weakly to fi^. 

Proof. In order to prove the first and second items, it is enough to use the preceding Lemma, 
Corollary 6 and Theorem 2 in |15j . plus the fact that 

hm NE[Af,{ijf] = -^B^-^ 
Af— >oo 2 — a 

. The third item is a consequence of the estimate 

C(m) < Cp™ 

with p = (^^i?") 2 and of Proposition 10 in [15J. The fact that /i^ is independent of the skewness 
parameter 9 is obvious since its moments only depend on the Cm for even m's, which are insensitive 
to the parameter 9. The fourth item is a consequence of Proposition 9 and Proposition 12 of |15j . 
The fifth one is a consequence of Theorem 1 of [15]. □ 

This lemma proves the first part of Theorem II. 81 In order to prove the second part we simply 
remark that we have already done so, since we have seen, in the proof of Lemma 13. H that fi^ 
converges and that its limit is the weak limit of E[/i^j^]. 



10 Appendix: Convergence to stable distributions for triangular 
arrays 

We begin here by recalling the notations for stable distributions, see for instance |llj . A real 
random variable Y has a stable distribution with exponent a G (0,2), a ^ 1, scale parameter 
(T > 0, skewness parameter /3 G [—1, 1], and shift parameter /x G M (in short Stablea{o', P, n) ) iff 
its characteristic function is given by: 

vro; 

E[ex.p{itY)] = exp [-cr''|t|"(l - i(3sign{t) tan(— )) + ifit] 

We will consider here only the case where a < 1. 

A complex random variable Y has an a-stable distribution with spectral representation (F, /i) if 
r is a finite measure on the unit circle S^, and is a complex number such that the characteristic 
function of Y is given by: 

E[exp(i(t,y))] =exp[- / s)|"(l - isi5n((t, s)) tan( — ))r(ds) + t)] 

We will need the constant 

inx , r(2-a)cos(^) 

ax = — 

X" 1 - a 

Throughout this section, we consider a sequence of i.i.d non negative random variables (X^.)fc>i and 
assume that their common distribution is in the domain of attraction of an a-stable distribution, 
with Q G (0, 1), i.e that the tail is regularly varying: 
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We introduce the normalizing constant a^v by: 

UN = mi{u,P[X >u]<^) (50) 

We consider a triangular array of real or complex numbers (Gjv,A;, 1 < ^ < N) and give sufficient 
conditions for the normalized sum: 

1 ^ 

to converge in distribution to a (real or complex) stable distribution. We will always assume that 
the triangular array is bounded, i.e that 

M := sup(|G7v,fc|,iV > 1,1 < k < N) < OD 

We begin with the case where the numbers Gjy^k are real. 

Theorem 10.1. Assume that the triangular array of real numbers {Gjy^k:^ > 1, 1 < ^ < N) is 
bounded. Furthermore assume that the empirical measure 



k=l 



converges weakly to a probability measure v on the real line. Then the distribution of the normalized 



sum 



Sn = z^Ylk=iGN,kXk converges to a Stablea{cr, P,0) distribution, with 



= ^ 1 \xrdu{x), 

J \x\°' sign{x)i^{dx) 



f \x\'^iy{dx) 

^N 

siaiemem snouia oj course oe unaersiooa as: ^N = l_ 

converges in distribution to zero 



Ifo'°' = 0, i.e iff = 6q, the above statement should of course be understood as: Sn = ^ Z]a~i GN,kXk 



Proof of theorem \10.1\ We begin with the particular case where the numbers GN,k are positive and 
bonded below. We assume that there exists an (5 > such that for any > 1 and 1 < k < N 

S < GN,k < M. (51) 

In this context we will be able to apply directly classical theorems to the array of non negative 
independent random variables 

UN,k = :^GN,kXk 
ON 

For instance, we could apply the theorem in section XVII. 7 of [6]. We rather choose to apply Theo- 
rem 8, chapter 5 of [7]. According to this last result, Theorem 1 1 . 1 1 will be proved in this restricted 
case if we can check the following three conditions. First the Uniform Asymptotic Negligibility 
(UAN) condition, for every e > 

lim max F(Un fc > e) = 0. (52) 
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Second we must check that: 

iUN,k<e)] - 0' (53) 



l<A;<Af 

and finaUy we must check that, for x > 



and that 



We first note that 



C 0"" 

hm P( max Um k < x) = expf ), (54) 



hm P( min UNk < x) = 1. (55) 

AT^oo \<k<N ' - ' ^ ^ 



'GN,k ~ (^)° 

which shows that ([5^ is thus a direct consequence of our assumption ([?T|) and of the fohowing 
lemma. 

Lemma 10.2. Let L be a slowly varying function and define ajsj as in (150)) ; 

aN = mi{u,P[\X\>u]<^) (56) 
Then , for any < a < b and any a < y < b 

^^ = ii(l + .(x.«)) (57) 



with 



hm sup e{x, N) = (58) 

N->oo a<y<b 



Proof of Lemma \10.B . Writing 

L{yaN) L{yaN) NL{aN) 1 



(yaiv)" L{aN) a% Ny^ 
this lemma is clearly a direct consequence of the classical fact: 



(59) 



hm = 1 (60) 

N^oo (oat)" 

and of the uniform convergence theorem for slowly varying functions (f2], Theorem 1.2.1), which 
asserts that the convergence 

is uniform for x's in a compact subset of (0, oo). □ 
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Next, in order to control the variance Var[UN,k^(UM k<'^)^ prove the validity of we must 
use Karamata's theorem, or more directly Theorem VIII. 9. 2 of [6jwhich shows that 

U„ = (62) 

Using this for Q = 1,2, we see that 

^-[f^^,a(^.,.<.)] - ^^'[yife]- (63) 



^ GN,k 
2 

Lemma (jl0.2p then shows that Var[Ui\f^k^(^if^ ^^^^] is of order ^, and thus that 



hm hm Yl ^^^pN,kliUj,,,<e)]=0 (64) 



l<k<N 



In order to complete the proof of Theorem llO.ll in the particular case where the numbers Gj\f^k are 
positive and bounded below, we now only have to check ()54p since ()55p is obvious, the variables 
U'N,k being non negative. For x > 



logP( max UN,k < x) = y^log[l 



1 ^GN,k^ 



l<k<N ' ^ ( ^f^N 

Using again Lemma (ll0.2p we see that 



C (T^ 

lim log P( max U N,k 1^ x) 



N^oo l<k<N X' 

since 



N 

f^oo iV ^ — ' J 



N 

1 

This checks the condition (I54p and finishes the proof in the particular case where the numbers 
G]y,k are positive and bounded below. Now it is easy to prove theorem 110.11 in full generality. It is 
enough to split the sum into the three independent summands 

^ N N 

Sn = }^ GN,kXk = \^ UN,k = Sti'^ — Sj^'^ + SN,e 

k=l k=l 



with 



TV 

= ^C/Ar,fcle<Gjv,fc 
k=l 
N 

= -'YUN,klGN,k<- 
k=l 

N 

SN,e = YUN,kl\GM,k\<<^ 
k=l 
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We now know that, if e and — e are not atoms of u, then S^'*^ (resp converges in distribution 

to a S'ta6/eQ,(cj^^, 1, 0) (resp Stablea{cr~ ^,1,0)) with 

/oo 
|x|"z/(da;) 



So that the sum S'^y'^ + 5*^'^ converges in distribution, when N tends to oo to a Stablea{o'a,e, l3a,e,0) 
with 

CaO'ae = / \x\"l'{dx) 
J\x\>e 

I\x\>e \x\°'sign{x)i^{dx) 



Pa 



lM>M"''(dx) 



It is clear that, since lirrie^o ^a,e = and that hmg_^o /5ci,e = Pa^ the distribution StablcaicTa^e^ Pa.e^O) 
converge to Stablcc^ ((Ta,/5o,0), when e tends to zero. Thus there exists a sequence eTv tending to 
zero such that the sum 5^'*^^ + 5"^'^^ converges in distribution to a Stablea{cra, (3a,0) variable. 
But 5'7v,ejv converges to zero in probabihty when N — > oo. Indeed, for any x > 0, 

1 ^ 

P(|%,^| > X) < VXfc > — ) 

aiv^ eN 

so that 

hm P(|%ejvl >x) = (65) 

These two facts show that = ^ Z]fcLi G^^k^k converge in distribution to a Stablea{cra, (3a,0) 
variable as announced in Theorem llU.il □ 

This result implies easily the following analogous result in the complex case. 

Theorem 10.3. Assume that the triangular array of complex numbers {Gj^^k,^ > 1,1 < ^ < N) 
is bounded. Furthermore assume that the empirical measure 

1 ^ 

k=l 

converges weakly to a probability measure v on the complex plane. Then Sn = 5^ 'l2k=i GN,kXk 
converges in distribution to a complex stable distribution with spectral representation (r,y,0) where 
Ti, is the measure on obtained as the image of the measure ■^\z\°'u{dz) on the complex plane 
by the map z ^ Again if u = 5q the above statement should be understood as: Sn converges in 
distribution to zero. 

Proof. For any fixed t G C, a direct application of Theorem IIU.II to the array of real numbers 
{{t,Gj\f,k)) shows that {t,SN) converges in distribution to a Stablea{cr{t), P(t),0) variable, where 
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^(*r = \{t,z)rdu{z) 

and 

^ J\{t,z)\°'sign{t,z)du{z) 
I\{t,z)\-di.{z) ■ 

As a consequence, we obtain that 

lim E[exp{i{t, Sn))] = exp - ip{t) tan(^))] 

Note that, by definition of T^,: 

a(tr(l-z/3(t)tan(^))= [ \{t, s)\^{l - isign{{t, s))tM^mu{ds) 

These two last facts prove that the distribution of Sn converges to a complex a-stable distribution 
with spectral representation (T^,0). 

□ 

In Section [5] we need a slight variation of Theorem I10.3[ We want to extend it to the case 
where the random variables Xk are truncated at a high enough level. More precisely, keeping the 
notations and hypothesis of Theorem 110.3^ we define, for any 6 > 0, the truncated variables 

We then consider the normalized sum 

TV 



-<5 



Theorem 10.4. Assume that the triangular array of complex numbers {GN,k,N > 1, 1 < < A^) 
is hounded. Furthermore assume that the empirical measure 



= ^ £ '^^iv,fc (66) 
fc=i 

converges weakly to a probability measure uon the complex plane . Then Sf^ converges in distribution 
to a complex stable distribution with spectral representation {r,^,0) where Fjy is the measure on 
obtained as the image of the measure ■^\z\'^i'{dz) on the complex plane by the map z — > Again 
if V = 5q the above statement should be understood as: Sn converges in distribution to zero. 

The proof of this variant is identical verbatim to the proof of Theorem 110.3^ we omit it. 

Finally we also need in Section [5] an information about the Fourier-Laplace transform of certain 
complex stable distributions. Consider a probability measure on C and define as above the 
measure T^. Let us denote by P'^ the complex StablcaiJ'u,^) distribution. 
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Theorem 10.5. Assume that the measure v is compactly supported in the closure of C . Then, 
for any t > 0: 

J e-^*^dP^(x) = exp(-r(l - a){it)" J x'^dv{x)) (67) 

Proof. This is a simple consequence of the analogous result for real Stableairr, 1) distributions. If 
X is a random variable with Stableaic, 1) distribution, and if 7 E C with 3^(7) > 0, then 

E(e-^^) = e"^^^^" (68) 

This result is classical when 7 is real positive (see Proposition 1.2.12 of [llj for instance). The 
statement (|68p is obtained by an easy analytic extension from the real case. 

Consider now a sequence of i.i.d.r.v {Xk)k>i, with common distribution Stablea{cr, 1). Further- 
more consider a bounded array of complex numbers (G7v,A:) £ C~, such that the empirical measure 
J! ^Gn k converges to z/ when N ^ 00. As above define the normalized sum 

N 



G 

Then, if ^N,k = it ^'^ , one has obviously 



k=l 



N 



E(e-^*^^) = n IE(exp(-7^,fcXfe)) 

k=l 

Noting that ^{'yN,k) > 0, it is then possible to use (j68|) : 

N 



E(e-*^-) = exp(--^5]7.V) 



Using the classical tail estimate for real Stablcaia, 1) distributions, when u tends to 00: 



P(X > u) 



one sees that oa? ~ Ca N a . 



Thus, we get the estimate 



jj'}2k=i^% k converges to J x'^dv{x). Using now the convergence theorem 110.31 we see 

that, 

f e-**^'ciP'^(x) = lim E(e-^*^'^) = exp(-- 3— — (if)° ( x'^dv{x)). 

J N^co CaCOS( — ) J 

Noting that 

„ /Tra, 1 — a 1 

Cn COS 



2 ' T{2-a) r(l-a) 



proves Theorem 110.51 □ 
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